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ON LATERAL CABLE OSCILLATIONS OF 
CABLE -CONNECTED SPACE STATIONS 

By Willard W. Anderson 
Langley Research Center 

SUMMARY 

The influence of lateral cable oscillations on a spinning cable-connected space sta- 
tion is defined in terms of steady-state and transient response to torques acting on the 
space station. The system is defined by using linear rigid-body equations and a dual 
formulation of the cable modes. Damping is included in the dual formulation and esti- 
mates of this damping are provided. The analysis allows a determination of the required 
depth of modeling for a given system in terms of the system parameters and the expected 
disturbances. Energy flow within the system is discussed by using a system bond graph. 

INTRODUCTION 

The equations of motion of a rotating cable — counterweight — space-station system 
in a gravitational field gradient are extremely complex and nonlinear. A complete solu- 
tion of these equations is therefore not feasible and thus the requirement exists that 
assumptions be made which reduce the mathematical problem to a level where solution 
is possible. 

This paper represents an attempt to define one aspect of the system, namely the 
lateral oscillation of the connecting cable and in particular the effect of this oscillation 
on the dynamic response of the space station and counterweight. Coupling between 
lateral and longitudinal oscillation is not included. The influence of the gravitational 
field and of the angular momentum of the system is not directly included. The variation 
of cable tension is neglected and only linear equations are used in defining the solutions. 
Low cable curvature arises from high-velocity lateral cable waves. Cable lengths are 
assumed to be less than 2 statute miles (3.2 km) and space-station and counterweight 
mass, geometry, and inertia are assumed to be of equal order of magnitude. These are 
important simplifications (refs. 1 and 2) which can influence the results, but they have 
been made in order to define clearly the basic character of the interaction between lateral 
cable modes and the space station and counterweight and in particular the effect of lateral 
cable forces on the attitude of the space station and counterweight. The analysis and 
description are presented in terms of easily calculable frequency parameters which refer 



to uncoupled-cable and massless-cable systems. The dual modal formulation is 
required, rather than the formulations of earlier papers. 

SYMBOLS 


A 

A,B,C 

A,B,C,D,E,F 

A 1’ A 2 


c 

E 


E' 


E" 

e 


F 1’ F 2 

f 


Sn 


* sfejje 

^n’^n’^n 


I 


Ii,I, 


1 ’ 2 


cable cross-sectional area (see eq. (8)) 

X-, Y-, and Z-axis rigid -body inertias 
steady-state parameters 
cable attachment lengths 
cable phase velocity 
cable complex modulus 
cable storage modulus 
cable loss modulus 
power bond effort 

base of natural logarithms (see, for example, eqs. (C6) and (D4)) 

Z -direction disturbance forces 

power bond flow 

cable mode spatial functions 

cable mode time functions 

cable area inertia 

system inertias 



system inertias (see eqs. (D2)) 


J 1> J 2 

k p k 2 

k 


system stiffnesses 
cable wave number 


L,N,P,Q steady-state parameters 

l cable length 

M c cable bending moment 

M Q? MpM 2 ,M x? M z disturbance torques 
m bond graph modulus 

m^n^ system masses 

n modal index 


P-pP 2 X -direction disturbance forces 

r cable damping parameter 

r c cable radius 

s cable coordinate 


T cable tension 


TpTg disturbance torques 

t time 


U 1’ U 2’ U 3’ U T system energies 
u_(t) unit step function 

o 
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u l’ u 2 


system translations 


V cable shear force 

v o’ v l’ v 2 ,v * system translations 
X,Y,Z system axes 

x, y,z system coordinates 

a x ,a-y,a z rigid -body attitudes 

ri’ y 2» y 3» y 4 steady -state phase angles 
e cable ratio parameter 

r] lateral cable displacement 

6 q, Op 63 system attitudes 

P n cable modal damping coefficient 

pA cable mass per unit length (p in this case denotes mass density) 

a cable tensile stress 

<pp<£>2><p3 system attitudes 

if/ lateral cable displacement 

spin frequency or velocity 
uncoupled system frequencies (see eqs. (D 2 )) 
o> forcing function frequency 

o» n cable mode frequencies 

4 


1 man 11 11 ■■ 



Subscript: 


max maximum 

Dots over symbols indicate derivatives with respect to time. Primes indicate 
derivatives with respect to distance. 

DESCRIPTION OF PHYSICAL SYSTEM 

The basic configuration of the physical system is shown in figure 1. The three 
portions of the system are the space station (mass = m^, inertia = I-jJ, the counter- 
weight (mass = m 2 , inertia = Ig^, and the cable (mass per unit length = pA, 
damping parameter = r, tension = T). The X,Y,Z axes system is also shown. A 
cable coordinate s is assumed to be approximately equal to y. The axes system 
rotates approximately with the physical system about the Z (spin) axis at the angular 
velocity System geometry is defined by the cable length l and the distances A^ 

and A 2 from space-station and counterweight centers of mass to cable attachment 
points. The space-station and counterweight inertias about the X- and Z-axis are 
assumed to be approximately equal, within the stability assumptions of appendix A. 

Crew-motion disturbances are defined by the X- and Z-direction forces Pj, Pg, 
Fj, and Fg and by the torques M^, Mg, Tj, and Tg. Space-station and counter- 
weight attitudes and translations are (pg, u p and u 2 * or the XY plane (fig. 2) 

and 0g> v i> and Vg for the YZ plane (fig. 3). Lateral cable displacements are 

)//(s,t) for the XY plane and T?(s,t) for the YZ plane. 

The system is further described in a later section in terms of effort and flow power 
variables by using bond graph techniques. 

DISCUSSION OF EQUATIONS OF MOTION 

If the effects of cable mass and cable damping are negligible, the equations required 
to define the motion of two rigid bodies (space station and counterweight) connected by a 
cable contain 12 variables (if positive cable tension is assumed and constraints are neg- 
lected). The 12 variables represent six degrees of freedom for two rigid bodies. These 
equations would be highly coupled and nonlinear. However, if the motion of the bodies 
relative to the rotating axes system of figure 1 is small and the system rotates slowly, 
the nonlinearity and coupling can be neglected. With this simplified system description 
the effects of cable mass and damping can be easily added. This constitutes the approach 
used in developing the system equations presented herein. 
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In appendix A the rigid-body Euler or low-frequency characteristics of the system 
are discussed and basic inertial requirements are placed on the system from stability 
and '’stiffness" considerations. The conclusions of appendix A are that the rigid-body 
constant -spin -velocity assumption is valid, that the Z-axis rigid -body coupling rate is 
small, and that the rigid-body X-axis attitude can be described with a stable second-order 
equation. 

Appendix B contains derivations of the system equations of motion where the cable 
lateral forces are considered in a general form. The derivation for the YZ plane 
begins with the rigid-body equation for X-axis attitude, the first of equations (A4). Dis- 
turbances and equations of motion for the Y-axis are not included since they do not 
directly influence the lateral cable motion. The equations of motion for the XY (spin) 
plane are shown to be identical to those for the YZ plane for high-frequency response, 
for systems where the flexible-body mode frequencies are well above the spin 
frequency f2 0 . 

Appendix C contains a derivation of the cable equations of motion and a dual formu- 
lation of these equations in modal form, including equations for predicting the damping 
parameters for the modes as functions of the cable geometry and material. The final 
linear equations for the YZ plane which constitute a linear mathematical description of 
the physical system are 


m i*i = F i + T iff 0 .') 
m 2*2 = F 2 " T 0^’^ 

, 

1 lh + TAjSj = M^TAjffPU) 

I 2 § 2 + TA 2»2 = M 2 + TA 2|f 

T?(0,t) = v 1 + A 1 0 1 

7i(l, t) = v 2 -A 2 0 2 ) 

v(l, t) - T?(0,t) = W 3 


(1) 


(2) 
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Equations (1) represent the space station and countweight; equations (2), the condi- 
tions at the cable attachment points (no cable shear being assumed); and equations (3), 
the cable itself, both as a partial differential equation and as an infinite series of ordinary 
differential equations (modal description). 

The following sections present response data for steady-state and transient distur- 
bances for the YZ plane only. These data are valid for the XY plane also, for the 
frequency range where lateral cable oscillations are important, but are not valid for the 
XY plane for disturbances of extremely low frequency. For example, the resonance at 
u> = in evidence in the following steady-state data for the YZ plane is not present 
in the XY plane; but the resonances at co = anc * so forth are present. 

Steady-State Response 

The dynamic equations of motion in the YZ plane, equations (1) to (3), have been 
solved for a harmonic input moment = M 0 sin cot acting on the space station. This 
disturbance was chosen for two reasons. First, it directly causes large space-station 
attitude response and thus would influence pointing accuracies of any experiment (tele- 
scope, tracking experiment, etc.); secondly, it yields representative vehicle gain plots 
for any attitude-control-system analysis. The solution is summarized in appendix D, and 
a digital program listing and explanation for calculation of response information are 
included in appendix E. The program calculates vehicle gain only but can be easily mod- 
ified to calculate phase in order to obtain more control -system information by using the 
phase equations of appendix D. Vehicle gain is given in figures 4 to 12 in terms of the 
four dimensionless amplitudes: v i/ v o’ anc * v 2/ v o where 
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Mp 

TAj 



as functions of the eight dimensionless variables <x>j£lp 
1 ^ 2 , A 2 jl, and I 2 /J 2 w ^ere 


Pp Ui/top m l/ m 2’ 


a 


1 " 



1/2 



c = 



"1 


_ 7 rc 


l 


_ 7rrc 
Pl 2TZ 

J 2 = m 2 A 2 2 


The first dimensionless variable is the ratio of the disturbance frequency cu to 
the uncoupled resonant frequency £2-^. Plots of steady-state response data indicate that 
the basic flexible -body attitude response of the space station resembles the response of 
a second-order system, with the resonant frequency £2^. The second variable is 
the damping coefficient for the first flexible-body cable mode. Appendix C contains a 
derivation of this coefficient, and a discussion is presented later in this section. The 
third dimensionless variable o>i/£2^ is the ratio of the fundamental uncoupled cable 
frequency to the uncoupled space-station resonant frequency £2^. The separation 

of these two frequencies in the frequency domain is a measure of the coupling between 
space -station attitude and lateral cable oscillation. Appendix C also contains a deriva- 
tion of coj. The fourth variable £2g/£2j is the ratio of the uncoupled counterweight 
resonant frequency £2g to the uncoupled space-station resonant frequency £2^. The 
four remaining variables are ratios of mass, inertia, and geometry. 

The interrelationship between these variables and the system response character- 
istics is quite extensive in terms of possible data generation. For example, 4 ( 3 ^) or 
approximately 8800 response plots would be required for only three values each for the 
last seven dimensionless variables. For this reason, only the main characteristics of 
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the system are presented and a program listing and explanation are included in appen- 
dix E for systems other than those presented. Also, to avoid symbol repetition, the 
values of the last seven variables will be given as a set of values where, for example, the 
set (0, 50, 1, 1, 1, 0.01, 0.333) indicates that 


Pi - o 



50 



1 


m^ 

m 2 


1 



.= 0.333 


Figure 4 contains space-station and counterweight response data plotted as func- 
tions of the forcing frequency ratio for the preceding set of variables repre- 

senting a symmetrical system of the indicated geometry and with lateral cable modes of 
very high frequency. Coupling between lateral cable oscillation and space-station and 
counterweight oscillation is therefore negligible. The space-station attitude response is 
seen to resemble the response of a single -degree -of -freedom system, with a resonant 
frequency and an additional resonance at the spin frequency £2 Q . The spin- 

frequency effect is not a true resonance but represents a rotation or tilting of the actual 
spin plane. This effect is a rigid -body characteristic of the system for YZ -plane oscil- 
lation. It occurs at the spin frequency £2 0 , where the input moment is in phase 

with the rotation, becoming equivalent to a moment fixed in inertial space and therefore 
continuously rotating the system. The linear mathematics used in generating these 
steady-state response plots are not valid at the resonance but are valid near the resonance 
if the resulting amplitudes are small. Counterweight steady-state response is relatively 
small, with maximum amplitudes near co = £2 0 and oj - = £2 2 . Figure 4(b) is a plot 
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of space-station and counterweight translations, which are equal for this case since cable 
mass is negligible » fijJ. The low-frequency rigid-body response can be expressed 
mathematically by the equations of appendix A and the relation 


v l| ~| V 2|~( Z + A l + A 2 )“x 


The high-frequency response indicates a resonance at a> = but for freqi ncies other 
than this, amplitudes are small. 

Figure 5 contains four space-station and counterweight attitude plots which can be 
used to demonstrate the influence of the dimensionless variables on the separation 
between spin frequency and uncoupled resonance frequency. The equation that relates 
the two frequencies is 


% 
. £2i 


h M 1 + m 2V fi 2 


J 2 l 


m 


iA«i 



1 




Ao 

1 + — 

1 + ( 


f- 

> 

l 

V 

\ 


1 h 

J 


( 4 ) 


Equation (4) is the result of an algebraic manipulation of the following equation for cable 
tension (where cable mass is neglected): 


m.m, 0 

T = —(l + A, + (5) 

m l + m 2 ' 1 ° 

Equation (4) illustrates the dependence of the spin ratio on the variables I 2 ^J 2 , ^2 A ’ 
m l/ m 2’ ^j/^l’ an< *> * 0 a l esser degree, The plots of figure 5 illustrate that 

changing these variables will affect the lower resonance, according to equation (4), with- 
out appreciably changing the basic high-frequency character of the system, when the 
fundamental cable frequency is high and the spin ratio is sufficiently small 

(0 0 « Figure 5(a) is a response plot for the system (0, 50, 1, 1, 1, 0.01, 0.0333). 

The plot illustrates that a decrease in the inertia ratio I 2 /J 2 resu ^ s i- n a decrease in 
spin ratio without affecting the response near co ~ Figure 5(b) is a response plot 
for the system (0, 50, 1, 1, 1, 0.1, 0.333). The increase in A 2 /Z from the system of 
figure 4 increases the spin ratio, but again the response near co ~ £2^ is not affected. 
Figure 5(c) is for the system (0, 50, 1, 4, 1, 0.1, 0.333). The decrease in spin ratio 
relative to the results of figure 5(b) is a result of the increase in the ratio m^m 2 . 
Figure 5(d) is for the system (0, 50, 2, 1, 1, 0.1, 0.333). The increase in the ratio 
^2/^1 * s reflected not only by the increase in but also by the additional 

resonance at co ~ f2 2 * This additonal resonance has a pronounced effect on counter- 
weight response but only a small effect on space-station response. Figure 5 indicates 
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that if the spin ratio £2 q /£2j calculated by equation (4) is small, the system response 
can be considered as the uncoupled superposition of flexible-body motion on rigid-body 
motion (appendix A) near u> ~ £2-^ ~ (when » £2^). 

Figure 6 is an attitude response plot for the system (0, 50, 0.5, 1, 1, 0.01, 0.0333). 

This system deviates from symmetry by having — = 0.5, this condition being reflected 

£ 2 ^ 

by the strong counterweight resonance at u> ~ £2g. The influence of this third resonance 
on space-station attitude is not appreciable. Figure 7 is an attitude response plot for the 
system (0, 50, 2, 1, 1, 0.01, 0.0333). The counterweight uncoupled resonance £2g * s 
twice the space-station uncoupled resonance £2-^. Again the effect of the additional 
resonance is not appreciable for the space-station response. 

The previous data are for systems where the cable fundamental mode frequency 
is well above £2 0 , £2-^, and This is a required condition if a cable is to be modeled 

as a massless string, with no degrees of freedom of its own. Figure 8 contains plots for 

o>i 

the system (0, 2, 1, 1, 1, 0.1, 0.333), where — = 2. The decrease in ci>]/£2j implies 

Oj ' 

the addition to the system of an infinite number of degrees of freedom and therefore, for 
an undamped system, an infinite number of cable resonances. However, when the funda- 
mental uncoupled cable frequency is sufficiently high, the coupling, although in existence, 
will not appreciably alter response, other than to superimpose the cable modes at nearly 
the uncoupled cable frequencies co n . This is very nearly the case for figures 8(a) and 
8(b), although two new resonant peaks now appear on either side of the resonance at 
u> ~ £2^. Space-station and counterweight translations are plotted in figures 8(c) and 8(d). 
The two resonances near u> ~ £2^ are clearly defined, as is the resonance associated 
with The response for the higher cable resonances is too low to appear in the fig- 

ures. The space-station and the counterweight translations exhibit a response at several 
frequencies near £2^ and greater than £2j, where one translation (vj) is zero and the 
other (V 2 ) is finite. This condition seems more descriptive than significant. 

Figures 9 and 10 contain the same response information as figure 8, but for the 
systems (0, 1, 1, 1, 1, 0.1, 0.333) and (0, 0.5, 1, 1, 1, 0.1, 0.333), or for systems with 

CO -I C0 1 

— =1 and — — = 0.5, respectively. The figures illustrate how cable -to -counterweight 

and cable-to-space station coupling increases as the fundamental uncoupled cable fre- 
quency decreases. This result reflects more influence of the cable mass on the system 
response. The point is made that the system is quite involved dynamically if the 
uncoupled cable frequencies u> n are near the uncoupled space-station and counterweight 
frequencies £2j and therefore, any simulation of such a system would require the 
inclusion of the effects of cable lateral motion when the ratio co-i /£2i is not large. 
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Cable damping has been assumed to be small for all previous response data. This 
assumption is valid if the variable is small. The following order -of -magnitude 
analysis shows that for systems where the cable length is long, the damping coefficient 
p i will be negligibly small. In appendix C the modal damping of equations (Cl 3) is 
expressed as 


rco r 


P n = 


2T 


Substituting equation (CIO) yields 



( 6 ) 


(7) 


For a cable of circular cross section, the area inertia is given as 



( 8 ) 


where r c is the cable radius and A is the cable cross-sectional area. Substitution 
of equation (8) into equation (7) yields 



E 


" /r 


(9) 


f P | | 

where a = — is the tensile stress in the cable. It is assumed a priori that E is of 
A 

the same order of magnitude as o , this assumption being valid for most metallic cables. 
Therefore, the order-of -magnitude relation exists that 


p i 



( 10 ) 


This equation clearly supports the assumption of small damping since cable radius r c 
is several orders of magnitude lower than cable length l. The data presented are thus 
justified. 

Figure 11 is included for the system (0.1, 0.5, 1, 1, 1, 0.1, 0.333), for which it was 
assumed that another mechanism was available to provide cable damping and that this 
cable damping was mathematically similar to the damping p n . This system is identical 
to that for figure 10, with the exception that p^ = 0.1. Rigid -body space-station and 
counterweight responses are identical near u> ~ £2 0 . The absence of flexible-body vehi- 
cle motion in this frequency region is thus substantiated. For the frequency range where 
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u> > S2 0 , response is increasingly altered as resonant amplitudes for the frequencies near 
(jo ~ £2^ become smaller, until, for the high cable modes, resonances no longer exist. 
This effect is caused by the increasingly greater damping in each higher cable mode 
(eq. (6)). 


Transient Response 

Equations (1), (2), and (3) have been digitally integrated for two transient distur- 
bance moments acting on the space station. The first moment Mj = 100u g (t) represents 
the disturbance caused by a crew member who moves relative to the artificial gravity 
field in the Z -direction and then remains at a new position. The second moment 
Mj = 100u s (t)cos(l.lt) represents the disturbance caused by a crew member who moves 
continuously between two positions after an initial motion similar to that of the first case. 
These disturbances act on a symmetrical vehicle defined by the following coefficients: 


m l = m 2 = ^0 slugs (10 216 kg) 

Ij_ = I 2 = 40 000 slug-ft 2 (54 232 kg-m2) 
Aj = A 2 = 10 feet (3.048 meters) 
l = 2500 feet (762 meters) 


pA = 0.015 slug/ft 
T = 10 000 pounds 


(0.718 kg/m) 
(44 482 N) 


c = 817 ft/sec 


(249 m/sec) 


P 1 = 


10" 6 


These coefficients are expressed in steady-state notation as (lO - ®, 0.649, 1, 1, 1, 0.004, 
0.57l). Figure 12 contains steady-state plots for attitude and translation of the space 
station and counterweight. These plots illustrate the close proximity of the station, 
counterweight, and cable lower resonances in the frequency domain. The general char- 
acter of these plots is similar to that of the plots previously described. 


Figure 13 contains time domain plots of space-station and counterweight response 
for the step moment = 100u s (t). These plots contain transient motion superimposed 
on the steady-state motion predicted by figure 12 for a zero frequency or constant input. 


Figure 13(a) illustrates the transient overshoot of space-station attitude by a factor 


of two times the steady-state displacement 0 Q 


TA X 


10“3 radian 


and then subsequent 
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oscillation at the frequency fi^ = 


= 1.58 rad/sec. The space station then 



begins to radiate energy through the cable to the counterweight, with accompanying loss 
in amplitude. 

Figure 13(b) is a plot of counterweight attitude. An initial delay in response 
At = £■= 3.06 seconds is required for any energy to propagate from the space station 
along the cable to the counterweight. The counterweight then begins to oscillate at the 

uncoupled frequency fig = = 1-58 rad/sec, with increasing attitude, until a complete 

exchange of available energy has taken place between the space station and the counter- 
weight. The energy exchange then reverses. This process would repeat until damping 
had absorbed the available energy and the system had reached the steady state. (For 
the low damping of the system in fig. 13, this absorption requires significant time and 
therefore is not a noticeable effect.) 

Figure 13(c) is a plot of space-station translation v^. The low-frequency oscilla- 
tion is that predicted by the first of equations (A4) 

•• o M x 
+ % ^ = X 

namely 

M„ 

Ox = — ^s-(l - cos fi„t) 

Afi 0 2 


where 

M x = 100u g (t) ft -lb (M x = 136u s (t) N-m) 


A - If + Io + 


m 1 m 2 


1 2 rrir. 4. xxi. 


T 


or 


m 2 


+ m 2 


m 


l m 2 A 1 + A 2 + l 


■ll + A, + A 2 ) 2 = 0.222 X 10 10 slug -ft 2 ( 0.301 X 10 10 kg-m 2 ) 

, \V2 

= 0.106 rad/sec 


> ( 11 ) 


v^ = -v*(l - cos fi Q t) 
v 2 = v*(l - cos fi Q t) 


( 12 ) 
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where 


l + A, + Ap Ti» 

v* = ± ^r= 0.005 foot (0.0015 meter) 

2 Afi 0 2 

The value of v* can also be found from figure 12; v* ~ 0.5v o = 0.005 foot (0.0015 meter). 

Equations (12) describe the low-frequency oscillation of figure 13(d) as well. The 
time lag At = 3.06 seconds for counterweight motion to begin and the low-amplitude 
high-frequency oscillations illustrate small flexible-body motion and coupling effects 
superimposed on the rigid-body motion. 

The first example dealt with a simplified crew motion of constant magnitude and 
therefore the transient response of the system could be interpreted easily. The second 
example deals with an oscillatory crew motion, which yields a more complicated super- 
position of transient and steady-state response. 

Figure 14 contains time domain plots of system response for the oscillatory input 
= 100u s (t)cos(l.lt). Figure 14(a) contains space-station attitude. This response 
indicates two oscillations, one at the space-station frequency = 1.58 rad/sec and 
one at the forcing frequency u> = 1.1 rad/sec. The oscillation at the space-station fre- 
quency appears to be similar to that shown in figure 13(a), while the oscillation at the 
forcing frequency maintains a constant magnitude. The steady-state magnitude as found 
from figure 12(a) is = 1.70 o = x 10“ 2 radian, which is approximately the oscilla- 
tion magnitude found for 100 < t < 120 seconds, where the steady state was reached 
temporarily for the system of figure 13. 

Figure 14(b) also indicates two oscillations, the steady-state oscillation at the 
forcing frequency having the lower magnitude #2 = 0.22 0 o = 0.22 x 10 -2 radian 
(fig. 12(a)), in evidence for 190 < t <210 seconds. Therefore it may be seen that the 
oscillation at ^ = 1-58 rad/sec is more dominant than in figure 14(a). The same 

energy transfer between the space station and the counterweight occurs until system 
damping absorbs all available energy and the response is entirely steady state. 

Figures 14(c) and 14(d) are quite similar if the basic energy lag time 
At = 3.06 seconds is neglected. The high-frequency oscillation occurs approximately at 
the forcing frequency. The low-frequency beat may indicate the proximity of the forcing 
frequency u> = 1.1 rad/sec and the lowest cable mode u>i = 0.649^ = 1.027 rad/sec. 

The steady-state oscillation amplitude can be found from figure 12 to be 

V 1 ~ v 2 ~ 0*lv o = 0.001 foot (0.0003 meter), which is approximately the average amplitude 

shown in figures 14(c) and 14(d). 
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This discussion of transient response was included in order to examine, for some 
forcing functions, the deviation of actual response from that predicted by the steady- 
state nondimensional plots previously presented. Since all types of disturbance inputs 
are theoretically possible and since the system is an infinite-degree-of-freedom system, 
no positive conclusion about transient response is possible. However, it appears from 
the data, including that presented, that orders of magnitude of maximum attitudes and 
translations can usually be estimated by approximating the disturbance as a harmonic 
function (or series of functions) and doubling the steady-state maximum response. Also 
the effect of energy flow between the space station and the counterweight will cause a 
transient oscillation of similar magnitude to appear when the steady-state prediction for 
that portion of the system is small. A second conclusion about the transient behavior is 
that the effect of the cable lag time At is not, in itself, critical for systems typical of 
those presently conceived. 


Bond Graph Representation 

A bond graph for the system described by equations (1) to (3) is included (fig. 15) 
to provide additional description and to serve as a graphical tool for the discussion of 
energy flow within the system. References (3) to (5) contain detailed discussion of the 
elements of a bond graph. A simplified discussion is also included in appendix F. 

Energy enters the system at the effort sources of bonds 1, 5, 13, and 16. These 
effort sources represent the space-station and counterweight forcing functions Mj, Fj, 
M 2 , and F 2 which arise from crew motion, aerodynamic drag, and so forth. This 
energy then flows from the space station or counterweight into the cable, causing the 
dynamic motion of the system and the interaction between the cable, space station, and 
counterweight. 

Figure 16 contains time plots of the energy of the space station, cable, and counter- 
weight and the total energy for the system of figure 13. The total energy, plotted in fig- 

TA 1 (20 O ) 2 

ure 16(a), has a maximum value of = 0.2 ft-lb (0.27 joule), which occurs 

initially and represents a loss of artificial gravitational field potential energy when the 
crew member changes position in the field at time t = 0. This energy enters the system 
at bond 1 and appears as space-station kinetic and potential energy (fig. 16(b)). It then 
flows through bond 4 to other elements of the system until at 100 < t < 120 seconds the 
counterweight energy (fig. 16(c)) is equal to the space-station energy. At this point the 

TAi0 o 2 

available space-station energy lenergy above the steady-state energy ^ = 0.05 ft-lb 

(0.068 joule)\ is zero and the counterweight begins to radiate its energy through the cable 
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to the space station, the energy exchange thus being reversed. It is interesting to note 
that the maximum energy contained in the cable (fig. 16(d)) is only 4 percent of the maxi- 
mum energy present in the system. Also the magnitude of the total system energy is 
extremely low. This result indicates that a probable low-power attitude control system 
could provide active damping for the space station (for crew-motion magnitudes dis- 
cussed herein). 


CONCLUDING REMARKS 


In the preceding analysis, equations for the dynamic system have been developed 
and solved. The results of this analysis allow certain concluding remarks concerning the 
behavior of the system, within the boundaries of all connected simplifying assumptions. 

First, the basic attitude response of the space station is that of an undamped second- 
order system, with a resonant frequency Coupled to this basic response are rigid- 

body characteristics, usually in evidence at frequencies lower than and cable lateral 
mode effects, usually in evidence at frequencies higher than S2j_. The entire system can 
be considered to be the superposition of uncoupled rigid-body, flexible -body, and "heavy 

^ 1 D 0 

string” motion when — « 1 and — — » 1 where — is the ratio of spin frequency 

Sij 

Ui 

to uncoupled space-station frequency and — is the ratio of fundamental cable mode 

frequency to uncoupled space-station frequency. However, the example system does not 
meet these conditions. Damping from mechanisms involved with the viscoelastic bending 
of the cable is extremely small because of the low dynamic curvature of the cable. The 
low cable curvature arises from the high velocity of cable lateral waves. 


The counterweight has a weak effect on the steady -state response of the space sta- 
tion. However, for any real situation where the transient response will usually dominate, 
the counterweight is seen to absorb energy from the space station and then radiate this 
energy through the cable back to the space station. Therefore, in this sense, the contri- 
bution of the counterweight is in no way negligible. 

Finally, the orders of magnitude of the maximum attitudes and translations of the 
space station and counterweight could be estimated, in the cases considered, by approxi- 
mating the disturbance in question as a harmonic function (or series of functions) and 
doubling the steady -state maximum response. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 22, 1968, 
125-19-03-10-23. 
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RIGID-BODY CHARACTERISTICS OF THE ROTATING SYSTEM 

This appendix presents a solution of Euler's equations for the system, where the 
cable length is large and the overall system inertias are constant (fig. 17). If the system 
is assumed to have only low-frequency disturbances (o> « fij) acting on it, it will behave 
as a rigid rotating body. If, in addition to the rigid-body motion, the total disturbance 
profile is assumed to cause only small flexible-body displacements relative to the space- 
craft dimensions, the rigid-body inertias will be essentially constant and the total space- 
craft motion will consist of the superposition of small linear flexible -body motion on 
rigid-body motion. 

Euler's equations are 


M x — Aflix + (C - B)aya z 

O = Bay + (A - C)a x a z > (Al) 

M z = C'a z + (B - Aj&x&y j 

where a x , a y, and a z refer to the rigid-body rates, a x and ay are small 
attitudes of the system about the x- and y-axis, and A, B, and C refer to the X-, Y-, 
and Z-axis rigid -body inertias. Torques M x and M z are assumed to be small and to 
act along the X- and Z-axis, respectively. 

These equations are simplified by the assumption that the cable is long and therefore 
Y-axis inertia B is small, and that A and C are approximately equal. Also the spin 
velocity a z is assumed to be nearly constant and of magnitude f2 0 . Direct substitution 
of these assumptions into equations (Al) is not possible, however, since it would violate 
the physical necessity that 


A + B > C 

A + C > B > (A2) 

B + C > A 

J 

The inequalities (A2) can be rewritten as 
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_1 < c ~ A < 1 
B 


(A + C > B) (A3) 


Equations (Al) are now rewritten with simplifying assumptions as 


a x + C A ^ 0 2q! x = 


a y = C ^o<*x 


& z = ^ n o“x«> 


(A4) 


C - A 

where all initial values other than a z = £2 0 are zero and the quantity — — — is 
bounded between 1 and -1. 

The stability of the system is now apparent (and well known), namely C > A. If 
maximum stability is assumed the condition that must be met is obviously 


C - A _ x 
B 

which reflects the system maximum stiffness under the disturbance torque M x . This 
condition is used in the derivation of the equations of the flexible system. 

The magnitude of the Y-axis rate as given by the second of equations (A4) is seen 
to depend on the X-axis displacement and the Z-axis rate. For long cable systems the 
displacement a x is small and therefore the Y-axis rate is also small relative to f2 0 . 
Also the Z-axis acceleration is seen by the third of equations (A4) to be second-order 
small since the quantity M z /A is small. This result substantiates the constant-spin- 
velocity assumption. 
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LINEAR SYSTEM EQUATIONS 

The linearized Euler equation for rigid-body motion in the YZ plane (the first of 
eqs. (A4)) for a condition of maximum YZ-plane stiffness is 

a x + Sl 0 2a x = ^ (Bl) 

A 


This equation can also be derived from conservation of angular momentum according 
to figure 17(a) where the system is shown as a rigid body and the body forces from the 
artificial gravity field acting on the space station and counterweight are equal in magnitude 
to the tension T and directed normal to the spin axis. The equation is 




T(Z +A 1+ A 2 ) 


Ol-v = 


Mx 

A 


(B2) 


where 


A = 


m l m 2 
mj + m£ 


(l + Aj + A 2 ) 


2 


Ij + I 2 « A 


(B3) 


Equation (5) defines the tension as 


m-. m 0 , . . o 

T = 1 — —(l + Ai + A 2 )f2 0 ^ 

m l + m 2 


Substitution of equations (B3) and (5) into equation (B2) will yield equation (Bl) and thus 
verify the uncoupled rigid-body assumptions of figure 17(a) and equation (B2). 

In a manner analogous to that for the derivation of equation (B2), equations for the 
complete system in the YZ plane can be derived. According to figure 3, conservation 
of lateral momentum requires that 
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while conservation of angular momentum requires that 



Also the constraints exist that 

>7(O,t) = v 1 +A 1 0 1 ^ 
V(l ,t) = Vg - A-2&2 ^ 

m, t) - T?(o,t) = i 


(B5) 


(B6) 


A comparison of figures 2 and 3 shows the similarity of the XY- and YZ-plane con- 
figurations, the exception being the direction of the body forces T (fig. 17) and thus the 
incorporation of a rigid-body oscillation at the spin frequency f2 0 . The response in the 
XY plane is therefore identical to that in the YZ plane for disturbances at frequencies 
much greater than the spin frequency £2 0 . 
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CABLE EQUATIONS AND MODAL DESCRIPTION 

The differential equation of motion for the cable is derived by considering conser- 
vation of lateral momentum for each incremental identified cable mass point. The cable 
is an infinite-degree-of-freedom system since there are an infinite number of mass 
points. An arbitrary mass point is shown in figure 18. The D'Alembert loading 
pA ds, the cable tension T, the cable shear force V, and the cable bending moment 

at 2 

M c are also shown in the figure. From lateral (^-direction) equilibrium, for small 77 , 


8 _ rp M .. SV^ds = pA ^5 ds 


“ rp 

as as as 


) 2 r 

at2 


(Cl) 


From rotational equilibrium (the rotary inertia being neglected), for small 77 , 


8M, 


as 


— ds = -V ds 


or 


ay _ 9 2 M c 
9s " a s 2 


(C2) 


Substitution of equation (C2) into equation (Cl), for constant tension T, yields 

82^ _ 9%c _ A 9272 

9 S 2 as 2 p at 2 


(C3) 


The cable moment M c is small but is included to introduce damping (energy dissipation) 
into the system. This moment is related to curvature by the equation 

M c = El ^ (C4) 

ds z 

If the material is assumed to be viscoelastic, the modulus is complex, or 


M c 




(C5) 
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where E" represents the viscous portion of the complex modulus (loss modulus), E' 
is the elastic portion of the modulus, I is the area moment of inertia for the cable cross 
section, and u> is the oscillation frequency for 77 , which is necessary for use of the 
complex modulus but is not a true material property. The effect of bending energy stor- 
age E’ will not be included hereafter since damping only is the desired second-order 
effect. 

The solution of equation (C3) with the cable moment described by equation (C5) is 
an undesirably complex way to introduce cable damping. Therefore, an alternative way, 
valid for small damping only, is presented. It is noted that this approximation deals 
mainly with the second term M c in equation (C3). 

A wave traveling in the cable could have the form 

?7(s,t) ~ ek^t - ks) (C6) 


where is the oscillation frequency and k is the wave number. Substituting equa- 
tions (C6) and (C5) into equation (C3) and neglecting E' required that 


k = ± 


-1 ± 1 + 


4eu> 


2\l/2 




2e/c 2 


1/2 


(C7) 


where 


c 



e 



The modulus of the complex quantity e/c% is small and can equal zero. There- 
fore, equation (C7) is not in a desirable form for computation. A more acceptable form, 
which considers e/c^ as a small number is 
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k = 




~\ 




k = 





(C8) 


These four approximate solutions are necessary to meet the four required boundary 
conditions of equation (C3) when the moment M c is defined by equation (C5). However, 
for small e the second two solutions result in a very local exponential correction of the 
equation of motion of the cable near the cable attachment point. The length of this region 
is given by the equation (ref. 6) 


h 


e 1/2 


Therefore, by assuming the condition l e /l « 1 is satisfied, only two boundary conditions 
are sufficient. The wave-number equation is therefore taken as 

k - ± |(x 4- (C9) 

where rco, defined by the equation 

ro> =E”I^= E”IpAa->2 (cl0) 

c 2 T 

is a pseudo -damping parameter. The choice of r is somewhat arbitrary since E" 
is a cable parameter which varies with cable tension (i.e. , strain level) T, cable mass 
per unit length pA, and frequency co. 

The definition given by equation (CIO) of the damping parameter r also is con- 
venient if a modal -equation approach is used. The reason is that equation (C3) assumes 
the form 


T + r —^L = pA — 7 : (Oil) 

9s z 8s Z 8t 8P 2 

which has a separable solution. (Substitution of equation (C6) into equation (Cll) will 
yield the wave -number equation (C9) and thus verify equation (Cll).) Solving equation 
(Cll) by using cable modes is standard and will be outlined only. 
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A separable solution for ?7(s,t) is assumed to be of the form 


%(s,t) = g n (s)h n (t) 


(Cl 2) 


Substitution of equation (Cl 2) into equation (Cll) yields the two separated equations 


^n + ^Pn^n^n + w n - 0 
§n + Sn = 0 


(C13) 


where 

^> n mode frequency 

p n mode damping coefficient, 

c undamped phase velocity, 

For the boundary condition at the cable extremities (s = 0 and s = l ) of zero lateral 
force or 


2T 

t y/2 

PA/ 


T g(0,t) = T ga, t) = 0 


the solution of the s -dependent function becomes 


g n (s) = cos 


ms 

l 


(C14) 


with the condition that 


The following lateral forces at the cable extremities (numerically designated according 
to fig. 15) are introduced: 
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e 8 = -T |2(0,t) 
e 9 = T |? C > t) 


(Cl 5) 


Equation (Cll) is then rewritten by using the 6-function 6(x) and becomes 


T^+r - PA ^ = -e 8 C(0) - e Q 6 (l) 


ds' 


ds^dt 


at* 


'9 l 


(C16) 


The solution is assumed to be expressible as a series of normal modes of the form 

OO 

V(s,t)= ^ h n (t)cos2|£ (C17) 

n=0 " 

Substitution of equation (C17) into equation (C16) and use of the orthogonality of the modes 
yields the following uncoupled equations for h n (t): 



e 8 + e 9 

pAl 


h n + 2 Pn w n h n + 


u) 2 h 
n n 


e 8 + (-l) n e 9 

pAl/2 


"N 

(n = 0) 


(n — 1, 2, . . .) 


7(0, t) = f g = h Q + ^ h n 
n=l 


} (C18) 


OO 

7?(Z,t) = fg = hQ + ^ (-l) n fl n 

n=l 

The preceding mathematical description of the cable is standard and normally ade- 
quate. However, the requirement exists in any computer solution of a system that vari- 
ables must be integrated and not differentiated with respect to time. In this specific case 
a dual formulation of the cable mathematics is thus required, where flow variables 
(velocities) are integrated and summed to yield effort variables (forces). The dual formu- 
lation of the equations can be found by an algebraic inversion of equations (C 18) in Laplace 
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notation, where all initial conditions are set equal to zero. References 7 and 8 contain 
more specific information on the dual formulation for elastodynamic problems. Ref- 
erence 9 discusses the inversion problem. The algebra of the inversion is omitted and 
the formulation is given as 

"N 

*0 = T( f 8 - f 9 ) ( n = °) 

hn + 2 Pn w nhn + «nM = ^[fg - (-1)%] (n = 1, 2, . . .) 

e 8 = ^0 + ^ ^n 
n=l 

e 9 = -hj - 7 (-l) n h* 

n=l J 

Equations (C19) can be rewritten, where cable slopes are functions of attachment- 
point velocities and positions, as 


> (C19) 


n=l 

fj(^t) = e 3 - f (-D n hr 

n=l 

h**+ 2p n o) n h**+ cv n 2 h**= | (0,t) + (-l) n+1 7ft,t)] 

e _ 7ft, t ) - V(0,t) 

3 l 

These equations are an equally valid description for XY-plane lateral cable motion where 
7 — xj/ and 9 — cp. 


> (C20) 


(n = 1, 2, . . .) 
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STEADY-STATE EQUATIONS 


The governing equations for the space station and counterweight, equations (B4), 
(B5), and (B6), are rewritten as two boundary -condition equations for the cable. If 
Fj = F 2 = M 2 = 0 these equations are 


where 



For a harmonic disturbance torque 


= M q sin cut 


(D3) 


a solution for cable lateral displacement rj which satisfies equations (Dl) as well as 
equation (C9) will be of the form 

? 7 (s,t) = aje^ ; 2 s S j[ n ( Cl ;t + kis) + bje ^ sin(o»t - kjs) 

+ a2e^2 s CO s(cjt + kis) + b2e~^2 s C os(cot - kis) (D4) 


where 


ki = Re(k) 
k2 = -Im(k) 
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The coefficients aj, a 2 , bj, and b 2 are found by substitution of equations (D3) 
and (D4) into equations (Dl), which yields the matrix equation 


1 + B 

1 - B 

-A 

A 


a l I 1 


r- "\ 

Ci/z 

A 

-A 

1 + B 

1 - B 


b 1 /l 


0 

1 - D + CF 

(1 + D + CF)E 

C - (1 - D)F 

(-C + (1 + D) F)E 

< 

a 2 1 1 

> = < 

0 ^ 

-C + (1 - D)F 

(C - (1 + D) F)E 

1 - D + CF 

(1 + D + CF)E 


b 2 /« 


0 

_ . 






v_ J 


L J 


where 

D 1 

A = Dikj = -y^k^ 
D 1 

B = Dik 2 = -J- k 2 l 
Do 

C = D 2 k x = ki l 

D 2 

D = D 2 k 2 = -± k 2 l 

E = e‘ 2k * 

F = tan kjZ 


Cl 


MqAj 



D 1 = - 



, h »i 2 ii\ 
J 1 u 2 J\) 



°2 




a> 2 J 2 y 
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The displacements of the space station and the counterweight can now be written as 
solutions of the preceding differential equations as follows: 


01 = 


_ M o/ ] 


Kl sin cot + - 1^(0, t) 


1 - 


or 




jsL 8s 

fi! 2 






1 - 


or 

f2r 


By introducing the notation 


v 2 \ A 2 


CO* 


■— ^(0 ,t) = L sin cot + N cos cot 
8s ’ 


^ (i,t) = P sin cot + Q cos cot 
8s 


the displacements can be rewritten in dimensionless form in terms of maximum amplitudes 
and phase angles. Therefore 


where 


SjTAj /GiTAi 


M r 


M r 


sin 


max 


(cot + yj 


faTA^ m \ 

\ ' max 


1 +■ 


TA] 

Ml 


L + 


'TA] 


N 




= tan" 
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where 



where 


and 


where 




L, N, P, and Q are defined as 


L = k 2 (a x - b x ) - ki(a 2 - b 2 ) 
N = k^aj - b {j + k 2 (a 2 - b 2 ) 


and where 
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P = sin k^Z 
+ cos k^Z 


-(ajki +a2k2)e k 2 Z - (bjkj +b2k2^e k2 ^J 
(aik 2 - a 2 ki)e k2Z - (bik 2 - b2ki)e ^ 

Q = cos k^Z ^a-jkj + a2k2^ e 2 - ^bjki + b 2 k 2 je k2 ^J 
+ sin kjZ ^4k 2 - a 2 k^ e k 2 Z + (b^ - b 2 ki) e k ^j 


There are eight basic independent dimensionless variables required to define the 
system maximum amplitudes. These variables are rco/T, col/ c, A ijl, 

A2/Z, wjtl 2, and I2/J2. The required program input must be in terms of the following 
eight dimensionless variables, written here as functions of the preceding eight equation 
variables: 


— - — 

_ 77 ro> c 
Pl 2 T col 

— = 77 — — 

f2j o)l 
_ n 2 o) 

m l _ J 1 hf n 2^ fp) a 2 l 

m 2 II J2V W / \ fi l/ l Aj 

h [uLf(2*?h _L 
12 \^l) \ W / l A 2 

A2 _ A2 

1 T 

h = h 

J2 J 2 

The required output are the four dimensionless maximum amplitudes 
(^lTAl/M 0 j max , ( v iT/M 0 ) max , (^2TA^/M 0 ) ma x) and ( v 2T/M 0 J max . 
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COMPUTER PROGRAM FOR STEADY-STATE SOLUTION 

The computer program for steady -state solution of the system equations for an 
input moment Mj acting on the space station is listed herein.* Eight input variables 
are required for solution of the problem. These variables, in terms of the equation 
variables defined in the analysis, are 


DP (1) 
DP (2) 

DP (3) 
DP (4) 
DP (5) 
DP (6) 
DP (7) 
DP (8) 


oo 



P 1 




^1 

m 2 

h 

l 2 

i 



The results are printed out as the following dimensionless maximum amplitudes: 

VAR (1) = 


/V A i* 


M. 


VAR (2) = 



*The program discussion and listing which appear in this appendix were supplied by 
Cohoon and Heasley, Inc., consulting engineers, Cambridge, Massachusetts, under con- 
tract to Langley Research Center. 
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VAR (3) = 



max 


VAR (4) = 



The program is arranged flexibly to permit solution of the problem for many values 
of any or all of the variables in a single run. Specified in the input for each variable are 
the number of equally spaced values it is to have and the smallest and largest of those 
values. The program will solve the equations for each of the given values of DP (1) by 
using the first value of the other input variables. It then indexes to the second value of 
DP (2) and repeats the solution for all the values of DP (1). This process is repeated 
until the complete set is finished. Note that although this is a very simple way to obtain 
results, it is also a very easy way to generate a tremendous amount of data. 

The following discussion and program listing should be sufficient to permit anyone 
familiar with the computer to use this program. 


Input 

The input variables are read in on eight cards, each of which has four numbers on 
it. The numbers on the cards are read with the format (12, E15.8, 14, E15.8). The first 
number on the card specifies the variable. The other numbers on the card specify, in 
sequence, the initial value of the variable, the total number of values of the variable which 
are to be considered, and the final value of the variable. The following is an example: 


1 

0.10000000E0 

11 

10.0000000E0 

DP (1) 

Initial value 

Number of values 

Final value 


Data must be fed in with the first numbers in ascending order. If this is not done, 
a statement is made and the program does not continue. 

Output 

The printout of the input variables and the results is in the following form with an 
E15.8 format: 
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DP (1) = 


DP (2) = 


DP (3) = 

DP (4) = 


DP (5) = 


DP (6) = 

DP (7) = 

VAR (1) = 

DP (8) = 

VAR (2) = 



VAR (3) = 


VAR (4) = 



A line is skipped whenever such a set is printed out. 

There are a number of checks in the program that have certain printouts associated 
with them as follows: 

(A) DATA CARDS ARE NOT IN ORDER. This statement is printed out if the input 
cards are not in the correct sequence. When this is printed out, the program stops. 

(B) DP (3, 4, 5, 6) CANNOT BE ZERO. This is printed out if the input DP (3), 

DP (4), DP (5), or DP (6) is zero or becomes exactly zero during the run. The pro- 
gram does not stop but rather continues to the next index. If the variable is still zero, 
the statement will be printed out again and the program will index again. 

(C) WITH DP (1) = 1 EQUATION BECOMES INFINITE. This is printed out if 
DP (1) is equal to 1. This condition would set some of the equations in the program 
equal to infinity. When this is printed out, the program indexes and continues. 

(D) WITH EV (7) = 1 EQUATION BECOMES INFINITE. This is printed out if 
EV (7) is equal to 1. This condition would set some of the equations in the program 
equal to infinity. When this is printed out, the program indexes and continues. 

(E) MATRIX IS NOT SOLVABLE IN THIS FORM. This statement refers to the 
solution of the set of simultaneous equations in the program and may occur if the program 
is asked to calculate a resonant amplitude with insufficient damping. If the divisor is 
smaller than a test value, this statement will be printed out. Since computation error is 
increased with division by small numbers, this is a check to inform the person who ana- 
lyzes the data that numbers in the solution matrix are getting small. The test value in 
the program is 0.00000001. When this statement is printed out, the program indexes 

and continues. 
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Program Listing 


PROGRAM CPH( INPUT , OUTPUT ,T APE5= I NPUT , TAPE6=0UTPUT ) 

PROGRAM FOR MAX. OEFL ECT IONS FOR NASA 1 FORTRAN 4 C+H 
DIMENSION RINI 181 , JACK I 8) , RF l N ( 8 ) , OR ( 8 ) ,DELT ( 8 ) 

DIMENSION DP (8) , EV(8) , GO( 4 , 5 ) , COF (4) , AT 1 4, 4) , VAR t 4 1 , ER ( 4 ) 
lill DO 301 1=1,8 

READ ( 5,302) JOE , R IN I ( I ) , J ACK ( I ) , RF IN(II 
332 FORMAT ( I 2 , El 5. 8, I4.E15.8) 

IF ( JOE-I ) 303,304,303 

303 WRITE (6,305) 

335 FORMAT (28H DATA CARDS ARE NOT IN ORDER) 

GO TO 300 

304 CONTINUE 
301 CONTINUE 

DO 306 1=1,8 

IF ( JACK ( I )-l)307»308»307 
308 DE L T{ I ) =0 . 

GO TO 306 

307 DR ( I ) = J ACK ( I ) - 1 

UELII I ) = ( R F I N ( I ) — RINI ( I) ) / OR ( I ) 

306 CUNTINUE 

DP(8)=RINI(8) 

JACK8= JACK ( 8 ) 

DO 319 J 8= 1 , J ACK8 
DP(7)=RINI (7) 

JACK7= J ACK ( 7 ) 

DO 317 J 7=1 , JACK7 
DP ( 6 ) =RI N I (6) 

J ACK6= JACK (6 ) 

DO 316 J6= 1 , JACK6 
DP(5 )=RINI (5) 

JACK 5= JACK ( 5 ) 

DO 315 J5= 1 , JACK5 
DP(4)=RINI (4) 

JACK4= JACK (4 ) 

DO 314 J4= 1 » JACK4 
DP (3) = RINI (3) 

JACK3= JACK! 3 ) 

DO 313 J3 = l » JACK3 
DP(2)=RINI (2) 

J ACK2 = J ACK ( 2 ) 

DO 312 J 2= 1 , J ACK2 
UP ( 1 )=RINI (1) 
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JACK1= JACK ID 
DO 3 L i J 1 = 1 » JACKI 
TEST=. OOOOOOOl 

C BEGINNING OF SOLUTION WITH ABOVE INPUT VARIABLES 

00 1 1=3,6 
IF (DPI I ) ) 1,2,1 
2 WRITE (6,4) 

4 FORMAT ( 42H DP(3,4,5,6) CANNOT BE ZERO LOOK AT PAGE 2) 

GO TO 100 

1 CONTINUE 
PI=3. 14159 

EV(1)=DP( 1 )*DP(2 )*2./DP(3) 

EV{ 2)=PI*DP( 1) /DPI 3) 

EV(3)=0P(6)*DPt 7 )/DP(4)**2 
E V I 4 ) =DP I 1 ) 

EV ( 5 ) =0P ( 8 )*DP(4)**4/(DP(5)*DP(6) ) 

E V ( 6 ) =DP ( 7 ) 

EV ( 7 ) = DP I 1 )/DP I 4 ) 

EV ( 8 ) =DP ( 8 ) 

FE= . 5*AT AN (EV(1)) 

COi=EV(2)*(l./(l.+EV(l)**2 )**. 25 ) *COS (FE) 

CJ2=EV(2>*( l./( l.+EVI 1 ) **2 )**.25)*SIN (FE) 

IF (l.-EV(4)**2) 5,6,5 
6 wRITEt 6,77) 

11 FO RM AT ( 34H WITH DP ( 1 ) = 1 EQ. BECOMES INFINITE) 

GO TO 100 

5 CONTINUE 

CQN=EV( 3 ) *( l./(l.-EV{4)**2)) 

DON=-EV( 3)*( l.+EVI 5)-( l./DPI 1 ) **2 ) *E VI 5 ) ) / ( l.-EV(4)**2) 
IFIl.-EVI 7) **2)20, 21, 20 
21 WR I T E ( 6,67) 

67 FORMAT ( 34H WITH EV(7)=1 EQ. BECOMES INFINITE) 

GO TO 100 
20 CONTINUE 

D0R=-DP(7)*( (l.+DP(8)-(l./EV(7) )**2*DP(8) ) / ( 1 .-E V ( 7) ** 2 ) ) 

A=0GN*C01 

B=D0N*CO2 

C=D0R*C01 

D= D0R*C02 

ELN=2. 71828183 

E=ELN**(-2.*C02) 

F = SI N (C01)/C0S (COl) 

GO ( 1 , 1 ) = 1 . +B 
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GO 12,1) =A 

GOI3.1 )=1.-D+C*F 

GO 14,1 )=-C+l i.-D)*F 

GO 1 1* 2 ) = l «-B 

GO ( 2 « 2 } =— A 

GO ( 3, 2 ) = ( l.+D+C*F)*E 

GU(4,2)=(C-l l.+D)*F)*E 

GO { 1 » 3 ) =- A 

GO I 2 » 3 ) = 1 • +3 

GO ( 3 » 3 ) =C- ( l.-D) *F 

G0(4,3)=1.-D+C*F 

GOIl,4)=A 

GO I2,4) = l.-B 

GO I 3 » 4 ) = I -C + 1 l.+D)*F)*E 

GO 1 4 * 4 ) = { l.+D+C*F)*E 

GO ( 1 » 5 ) =CON 

G0(2»5)=0.0 

G0(3,5)=0.0 

GO (4,5 )=0.0 

BEGINNING SOLUTION OF MATRIX 
AT I 1 » 1 ) =GO ( 1 » I ) 

IF { AB S I ATI 1, 1 ) )-TEST)40,40, 16 
40 WRITE! 6, 11 ) 

11 FORMAT l 36H MATRIX IS NOT SOLVABLE IN THIS FORM) 
GO TO 100 

16 CONTINUE 

AT I 2 , 1 ) =G0 (2,1) 

AT 13 » 1 ) =GQ (3*1) 

AT { 4 * l ) =G0 (4,1) 

HI 2=G0 ( 1 » 2 ) / AT I 1,1) 

AT (2 , 2 ) =GQ (2,2)— AT{2,1)*H12 
IFIABS ( ATC2,2) )-TEST ) 40 , 40, 17 

17 CONTINUE 

ATI 3,2>=G0(3,2)-AT(3,l)*H12 
ATI4,2)=G0I4,2)-AT<4, 1 )*H12 
H 1 3= GO I 1 , 3 ) / AT ( 1,1) 

H23=(G0(2,3)-AT( 2, 1 ) * HI 3 ) / AT { 2 , 2 ) 

ATI 3,3)=G0(3,3)-AT(3, 1 ) *H1 3-AT I 3 ,2 ) *H23 
IFIABS I AT (3, 3) )-TEST) 40, 40,18 

18 CONTINUE 

AT I4»3)=GQ(4,3)-AT14» 1) *H1 3- AT I4,2)*H23 
H 14 = GU I 1 , 4 ) / AT I 1,1) 

H24=IG0I2,4)— ATI 2, 1 ) *H 14 ) / AT I 2 , 2 ) 
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H34=(G0( 3,4)— AT { 3, I ) *H L4-AT (3»2)*H24)/AT(3,3) 

AT (4,4) =GQ (4,4) -AT (4»1)*H14— AT(4»2)*H24-AT(4»3)*H34 
1FIABS ( AT(4,4) )-TEST)40,40, 19 
19 CONTINUE 

CRl=GO( 1, 5)/AT( 1,1) 

CR2=(G0(2,5)-AT(2,1)*CR1)/AT(2,2) 

CR3= ( GO t3»5)-AT(3»L)*CRl-ATi 3, 2 I *CR2 ) /AT (3, 3) 

CR4= ( GO (4,5) — AT (4,1) *C Ri-AT (4» 2)*CR2— AT (4, 3) *CR3 )/AT(4»4) 
GO ( 4 » 5 ) =CR4 

GO ( 3 , 5 ) =CR 3— H34*GO (4,5) 

GO ( 2 , 5 ) =CR2-H24*GQC 4*5) — H2 3*GO ( 3 » 5 ) 

GO( L ,5 )=CR1-H14*G0(4, 5 )-H13*G3 ( 3 , 5 ) -H 1 2*GO ( 2, 5) 

C ENO OF MATRIX SOLUTION 

C0F(1)=C02*( GO ( i»5)-GO(2»5))-CQl*(GO(3»5)-G0(4»5)) 
CQF(2)=C01*( GO ( 1»5)-G0(2»5))+C02*( GO (3,5) -GO (4,5)) 

X 1 = ( GO ( 1 , 5)*C01+G0(3, 5 >*C02) *ELN**C02 
X2=(-GO( 2,5)*C0l-G0(4,5)*C02)*ELN**(-C02) 

X3=( GU( 1, 5)*C02-G0( 3, 5) *C01) *ELN**C02 
X4=(G0( 2,5)*CO2-G0(4,5)*C01)*ELN**(-C02) 

COF ( 3 ) =S I N (COl) *(-Xl+X2)+COS ( CO 1 ) * ( X3-X4 ) 

CUF ( 4 ) =COS (C01)*(X1+X2)+SIN ( CO 1 ) * ( X3+X4 > 

CL 1=1. -DP( 1 )**2 

CL2= ( 1 . +COF ( 1 ) )**2 

CL 3= ( COF ( 2 ) ) **2 

VAR ( l ) = ( SORT ( CL2+CL3 ) ) /CL 1 

VAR(2)=(EV(5)/DP( 1)**2)*SQRT ( (COF(l) )**2+CL3) 
CUL=1.-EV(7)**2 
OL 1= ( C OF (3) ) **2 
OL2= (COF (4) ) 

VAR ( 3 ) = ( SORT ( DL1+DL2 ) I/COL 
DL 3= ( ( DP ( 7 ) / EV ( 3 ) ) *COF (3) ) **2 
DL4= ( { DP(7 ) / E V( 3 ) )*C0F(4) )**2 
VAR(4)=DP(8)*( 1 . / E V ( 7 ) )**2*SQRT (0L3+DL4) 

WRITE (6,600) DP( 1 ) »DP(2) ,DP(3) 

600 FORMAT ( 7H DP ( 1 ) = , E I 5 . 8 , 7H OP ( 2 ) = , E 1 5 . 8 , 7H OP ( 3 ) = , E 15 . 8 ) 
WRITE (6,601) 0P(4),DP(5),DP(6) 

601 F 0 RMAT ( 7H OP ( 4 ) = , E 1 5 . 8 , 7H DP ( 5 ) = , E 1 5 . 8 , 7H OP ( 6 ) = , E 1 5 . 8 ) 
WRITE (6,602) DP(7),DP(8) 

602 FORMAT ( 7H DP ( 7 ) = , E 15 . 8 , 7H DP ( 8 ) = , E I 5 . 8 ) 

WRITE (6,603) VAR( 1 ) , VAR( 2 ) 

603 FORMAT ( 10H VAR ( 1 )= , E 1 5 . 8 ,9H VAR ( 2 ) = , E 1 5 . 8 ) 

WRITE (6,604) V AR ( 3 ) , V AR ( 4 ) 

604 FORMAT( 10H VAR ( 3 ) = , E 1 5. 8 , 9H VAR ( 4 ) = , E 1 5 . 8/ ) 
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100 CONTINUE 

OPI I) = DPm+DELTI 1) 

31 1 CONTINUE 

DP (2)=0P(2)+DELT{2) 

312 CONTINUE 
DP(3)=0P(3»+DELT< 3) 

313 CONTINUE 
0P(4)=DPI4)+DELT(4) 

314 CONTINUE 
0P(5)=0PC5)+DELT( 5i 

315 CONTINUE 

OP (6) = 0P(6) +DEL T ( 6 ) 

316 CONTINUE 
DP(7)=QP(7)+DELT(7) 

317 CONTINUE 

OP (8J=0P(8) +OEL T ( 8 ) 
319 CONTINUE 

GO TO 1111 
300 STOP 
END 
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DISCUSSION OF SIMPLIFIED BOND GRAPH NOTATION 

This short discussion is intended to define the symbols used in figure 15 only and 
is in no way a definition of the complete notation of bond graphs. 

The power bond is defined with the symbol where the product of the signal 
pair e and f describes the power passing through the bond in the direction of the half 
arrow. For the system of figure 15 the effort e and the flow f are either the pair 
torque and angular speed or the pair force and linear speed. The signal pairs are not 
shown in figure 15 but are implied. A power bond number is assigned such that the effort 
associated with bond 1 is ej and the flow fj. 

The effort sources (e.g., E-^ of figure 15 denote that effort (e.g., ej) is a speci- 
fied function of time. The symbol 0 denotes a point of common effort; the symbol 1, a 
point of common flow. The storage element for potential energy has the symbol -*■ C; 
that for kinetic energy, the symbol -*■ I; that for energy loss, the symbol -*■ R. The final 

ci b 

element of figure 15 is the transformer - TF -, which implies the relations 

e a = me b 



where m is the transformer modulus. 
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(c) High mass ratio, m^m 2 . 
Figure 5.- Continued. 
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(c) Space-station translation. 
Figure 8.- Continued. 




(d) Counterweight translation. 
Figure 8.- Concluded. 
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(c) Space-station translation. 
Figure 9.- Continued. 
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(b) Counterweight attitude. 
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(c) Space-station translation. 
Figure 13.- Continued. 
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(d) Counterweight translation. 
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(a) Space-station attitude. 

Figure 14.- Transient response to oscillatory step input torque. 
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Figure 14.- Continued. 
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(c) Space-station translation. 
Figure 14.- Continued. 
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(d) Counterweight translation, 
Figure 14.- Concluded. 
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Figure 15.- Bond graph representation of system. (See appendix F for discussion of notation.) 





(a) Total energy. 

Figure 16.- Energy flow for constant step input torque. 
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(b) Space- station energy. 

50 Figure 16.- Continued. 
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Figure 18.- Incremental cable mass equilibrium. 
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